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Abstract. We have obtained new BVR-b&nd images of 
the gravitationally lensed QSO pair HE 2149-2745 A,B 
(QSO redshift z e = 2.033; angular separation 9ab — 
1'.'71), under 1" seeing conditions. Subtraction of a scaled 
point-spread-function reveals an extended object of an- 
gular size comparable to 6>ab, and R magnitude of 20.4 
lying between and equidistant to the QSOs. The geome- 
try of the system leads us to interpret this object as the 
lensing agent responsible for the double QSO images. A 
fit to the galaxy surface brightness profile indicates an el- 
liptical galaxy, which is further supported by the absence 
of Mg II in absorption. The lensing galaxy is not detected 
in B and V . We discuss possible lens models for this ge- 
ometry constrained by our observations and consequences 
for the mass and redshift of the lens. 
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Gravitational lensing of distant extragalactic sources (e.g., 
QSOs) by intervening large masses at low to moderate red- 
shift (e.g., galaxies) has become a powerful tool of mod- 
ern cosmology. The multiple images of lensed QSOs are 
expected to be separated by a few arc seconds in the sky, 
implying transverse line-of-sight (LOS) separations of up 
to few tens of kpc. The detection and study of the lensing 
agent provide an independent way to (i) probe the dark 
mass involved in the lensing potential (e.g., Natarajan et 
al. 1995); (ii) put constraints on the the Hubble parame- 



ter Hq in case of QSO variability (Refsdal 1964). However, 
the position of the lensing agent, normally much fainter 
than the lensed source, is fairly difficult to determine. To 
date, about 20 lens redshifts have been determined out of 
twice as many known gravitational lens events. 

In this letter we present new BVR images of the grav- 
itationally lensed QSO pair HE 2149-2745 A,B aimed at 
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detecting the galaxy responsible for the double image. The 
QSO redshift is z = 2.033 and the angular separation be- 
tween the QSO images is 6*ab = 1''71. As the QSO pair was 
discovered, the spectr oscopic and i?-band observations by 
Wisotzki et al. ( 1996 ) already revealed its lensing nature, 
but the detection of the lensing galaxy remained uncer- 
tain. Our new images confirm the presence of the lensing 
galaxy. The sampling is well below the small separation 
between the two QSO images, so an estimate of galaxy 
parameters is possible. They are discussed in the frame- 
work of possible lens scenarios. 

2. Observations and image analysis 

We observed HE 2149-2745 on August 8th, 1997, with 
SUSI on the 3.5m ESO NTT. SUSI uses a TEK CCD with 
1024 2 24 fim pixels corresponding to 0'.'13 in a 2'.2 x 2'2 
sky field of view. Despite a rather variable seeing several 
exposures could be taken under ~ 1" seeing conditions. 
This was lucky given the bad weather conditions just days 
before the observations. The night was not photometric. 

We obtained science exposures in the three bands B (5 
exposures a 600 s), V ( 5 a 400 s), and R (9 a 300 s). We 
also obtained short exposures of the photometric standard 
star field PG2213-006 (Landolt |i~992] ) to check for possible 
spatial variations of the point-spread function (PSF) and 
to normalize our magnitudes to the standard system. After 
bias-subtraction the images were corrected for detector 
pixel-to-pixel variations using combined twilight sky flat- 
fields. The effective width of a point source in the science 
frames, as measured for the only well-exposed star in the 
SU SI fie ld (hereafter "star 3"; cf. Fig. 1 in Wisotzki et 
al. 1996 ), range from 1'.'3 to L"5 in the B-band frames, 
from r/3 to l'/4 (V), and from l'.'O to 1'.'4 (R). Due to 
this quality dispersion we decided to analyse each image 
separately. 

All BVR exposures clearly resolve the QSO compo- 
nents (hereafter A and B, respectively). We used the 
DAOPHOT II package to derive the PSF and to obtain as- 
trometric and photometric parameters in the frames con- 
taining the QSOs. The PSF was modeled with a (3 — 2.5 
Moffat function whose parameters resulted from the fit to 
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star 3, allowing for empirical departures from the analyti- 
cal form. A similar one-star fit and subsequent subtraction 
of the scaled PSF in the more crowded standard-star field 
demonstrated the spatial stability of the PSF shape, thus 
validating the use of just one star in its calculation. 

3. Results 

3.1. R-band images: the lensing galaxy 

Fig. (left) shows an i?-band single exposure of HE 
2149-2745 A,B with FWHM = l'.'O. The subtraction of 
the scaled PSF at the position of both QSO images clearly 
leaves a diffuse, ellipsoidal object lying between and more 
or less equidistant to A and B, as is shown in the right- 
hand panel. The simultaneous scaling of the PSF overfit- 
ted the fluxes at the A and B peaks, thus leaving small re- 
gions with relative negative residuals around both QSOs. 
This is expected given the contribution from the underly- 
ing object outside a few seeing radii; consequently, small 
corrections of roughly 0.01 (A) and 0.02 (B) magnitudes 
were made to recover the lacking flux. The mean flux val- 
ues in the central region deviate 15a from the background 
sky level. The stability of the PSF leads us to conclude 
that we have undoubtedly detected a galaxy in the LOSs 
to HE 2149-2745 A and B, most probably the gravita- 
tional lens. 

Simple aperture photometry yields Rg — 20.35 ± 0.20 
for the galaxy. The error bar reflects the tolerance range al- 
lowed by the fine-tuning described in the above paragraph, 
and includes the zero-point uncertainty. For the sum of 
both QSO components we find Ra+b — 16.30 ± 0.04. 
These magnitudes are based on the photometric results 
by Wisotzki et al. (1996) for star 3 (assuming it has not 
varied), and consider a color-term correction of 0.05 which 
reduces them to the Landolt (1992) system^. 

The relatively deep _R-band images reveal many red 
non-stellar objects in the SUSI field of view, but no ev- 
ident overdensity is observed close to the lensing galaxy 
to a limiting magnitude of R ~ 23. However, we cannot 
exclude the possibility that the galaxy belongs to a cluster 
since the small field of view implies transverse distances 
still compatible with cluster sizes. In particular, two of the 
objects, not much fainter than 20.4, lie within 12" of the 
galaxy position. It remains therefore unclear whether the 
galaxy belongs to a cluster or not. 

Table g also gives the calibrated magnitudes of both 
QSO components. As the V magnitude of star 3 is not 
available we assumed V — R = 0.46 ± 0.06 using a selected 
sample of stars with similar £? — R colors. We note that the 
magnitude difference between A and B is consistent with 
a flux ratio A to B of / = 4.3 for all three band filters. 

Fig. U shows isophotal contours of the lensing galaxy. 
The isophotes are separated by 0.2 magnitudes, and the 
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Fig. 1. Left: i?-band 300 s exposure of HE 2149-2745 A,B. 
North is on the top, east to the right. Right: Same as in 
left-hand panel but after scaling and subtracting the PSF 
at the positions indicated with the crosses. Two exposures 
with similar effective seeing of l'.'O have been combined 
(600 s). 

Table 1. Photometry and astrometry of HE 2149-2745 
A,B and the lensing galaxy. 





A 


B 


m,A — me 


Galaxy 


X 


0.00 ±0.03 


0.90 ± 0.03 




0.50 ±0.08 


y 


0.00 ±0.03 


-1.45 ±0.04 




-0.73 ± 0.07 


B 


17.13 ±0.03 


18.70 ± 0.03 


1.57 ±0.04 


> 22. 2 a 


V 


16.82 ±0.06 


18.39 ±0.06 


1.57 ±0.08 


£ 21.7 a 


R 


16.61 ± 0.02 


18.20 ± 0.02 


1.59 ±0.03 


20.35 ±0.20 



1 The B magnitude for star 3 given in Wisotzki et al. (1996) 
was erroneous and must be corrected to B — 16.14 



Approximate 3<r magnitude limits. 



lowest contour level is approximately at 7a above the 
background count. Note the distorted region near QSO 
A, probably due to photon-shot noise near the A peak 
(see below). The positions of A and B are marked with 
filled squares, since their errors are too small to appear 
at this scale. The angular separation between A and B 
is 0ab = 1'.'71 ± 0'.'07, and the position angle is PAab = 
29.5° ± 1.0°. The position of the galaxy was derived using 
the isophotes at brighter levels than half the peak inten- 
sity: x = 0'.'50 ± 0'.'08, and y = -0'.'73 ± 0'.'07 relative to 
A. There is a small misalignment between this position 
and the QSO images, but alignment is possible within the 
measurement uncertainties. Its orientation coincides sur- 
prisingly well with the line joining A and B (AB), with an 
ellipticity of e = (a — b)/(a + b) ~ 0.33. The galaxy center 
lies nearly equidistant to both QSOs, with 9ag = CK'88 
and 6» BG = 0'.'83. 

The good spatial sampling of our data allows an es- 
timate of the galaxy parameters. A 180° section aligned 
with the south semimajor axis was used to fit isophotal 
ellipses to the flux values. Fig. || shows the averaged sur- 
face brightness profile (error bars). A de Vaucouleurs law 
and an exponential disk were fitted to the profile within 
0'.'5 < r < 1'.'7, thus including only values at the > 3a 
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Fig. 2. Isophote contours of the lensing galaxy in intervals 
of 0.2 magnitudes. The lowest contour level is at la above 
the sky background. 



level but not significantly affected by the seeing smooth- 
ing. The final fit models result from re-constructing the 
isophotal ellipses with the fitted profiles and smoothing 
them with the PSF. We see that the data are better mod- 
eled by a r 1 / 4 law (xt = 0.3) than by an exponential disk 
(x„ = 2.2), especially at the core. Furthermore, images 
taken under better seeing conditions should accentuate the 
concave shape of the observed profile. We take this result, 
though cautiously, as evidence for an elliptical galaxy. 

3.2. B and V band images 

The galaxy is not detected in the PSF-subtracted B and 
V band images. However, the region surrounding QSO 
image A is slightly overfitted at both sides of AB. Such 
residuals are not observed for B. This "symmetric" overfit, 
though not significant 2a in B and only marginal in 
V), suggests A might be composed of two or more fainter 
point sources lyi ng on AB. We inves tigate d this possibility 
(cf., Bade et al. |l997t Burud et al. |1998| ) in the S-band, 
as one expects here less contamination by a hypothetical 
foreground galaxy (in spite of the fact that the PSF is 
broader in the blue). However, attempts to re-fit A with 
two sources of nearly half its intensity failed at recovering 
the background level. It is difficult to establish on the basis 
of the present data what causes the low-quality fits in the 
B images. From an observational point of view, splitting 
of QSO image A in very close sub-components cannot be 
ruled out, but an underlying object could also contribute 
to slightly distort the QSO images. An explanation of this 
must await better quality data. 

To put upper limits on the galaxy B and V brightness a 
variance frame was created which considers photon statis- 
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Fig. 3. Surface brightness profile of the lensing galaxy. 
Solid line: r 1 / 4 law fitted to the data within CK/5 < r < 
l'/7, and convolved with the PSF. Dashed line: Same for 
an exponential disk. The fitting range is indicated by the 
vertical dotted lines. 

tics (dominated by the QSO fluxes), readout noise, and 
the uncertainties introduced by background subtraction. 
Integration of the variances in a 3'.'5 radius aperture yields 
detection limits for B and V. In addition, we consider the 
fraction of possible galaxy light hidden in the QSO see- 
ing disks and hence subtracted with the PSF. Analysis of 
R images with similar seeing as in B and V shows that 
this fraction can be as large as 0.7 (B) and 0.5 (V). The 
estimated (3a) detection limits are listed in Table [l| 

4. Discussion 

4-1. Redshift of the galaxy 

The redshift of the lensing galaxy can in principle be con- 
strained by its color if we know its morphology. Besides the 
surface-profile fits, further evidence that we are observing 
an elliptical galaxy in front of HE 2149—2745 comes from 
the optical spectra of QSO A: no Mgn absorption system 
is detected down to an (observed) equivalent width limit 
of ~ 0.3 A (3a). Given the small impact parameter (~ 7 
h^kpc at z — 1), this result almost excludes a disk-like 
galaxy in the foreground at z > 0.35. 

The 3cr lower limit derived for the galaxy V magnitude 
(Table 0) puts a lower limit of V - R ^ 1.4. This is the 
V — R color an E galaxy would have at a redshift of ~ 0.3 
or larger (using spectral energy distributions observed at 
z = from Coleman, Wu & Weedman 1980| ). A bound 
consistent with this redshift is obtained from the B — R 



color (Bressan, Chiosi & Fagotto 1994). The implied abso- 
lute luminosity for zq — 0.3 is Mr = —21.5 (considering 
if -corrections, Hq — 50 km s^ 1 Mpc^ 1 , and go = 0.5), 
i.e., very close to M* . An upper limit for zq is difficult to 
establish, but for the observed Rq the expected luminos- 
ity beyond z ~ 0.5 becomes too large to be real. We thus 
arrive to 0.3 zq ^ 0.5. 
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4-2. Lens models 

Because of the symmetry of the mass distribution ex- 
pected for any regular galaxy, the deflection angle and 
the amplification should be very similar for the two im- 
ages of HE 2149—2745. This is also true if we include an 
external shear. To explain the flux ratio of 4.3, the depen- 
dence of the amplification on the positions has to be very 
strong. This can be achieved if the images are located near 
a critical curve, implying high amplifications. 

Given the required sensitivity of the models for small 
changes of the positions and the small number of con- 
straints, a maximum likelihood model fitting is not ap- 
propriate for this system; instead, we use an analytical 
approach to find the possible model parameters consider- 
ing the measurement uncertainties. 

We use a singular isothermal elliptical mass distribu- 
tion (SIEMD) as given by Kassiola & Kovner Ql993p . As 
can be seen from Fig. |^, the images are almost exactly 
located on the major axis of the galaxy. To simplify the 
calculations, we use the line AB as the major axis and 
project the center of the galaxy onto this line. We further 
include an external shear 7, whose source has to be lo- 
cated on the major or minor axis to be in agreement with 
the observed image positions. As observational parame- 
ters, we use the ratio of distances of the images from the 
center of the galaxy f x = #bg /$ag (nearly unity) and the 
amplification ratio / = 4.3. In addition, we force the two 
images to have different parity, which is a necessary con- 
dition to exclude the existence of more than two images. 
Even non-singular models (PIEMD) rule out the possible 
splitting of A in the radial direction. 

On the main axis, the lens equation and the amplifi- 
cation for the SIEMD model with external shear read 



(l-j)x-a 
(1-7) fl + 7 



1 



e 2 2V~e • 
arctan sign x 



2^~e I < 

a (1 + e) 



Ml 



Since a degeneracy in the models prevents the inde- 
pendent determination of e and 7, we use the two above 
equations to define the parameter E, 



E = 



1 



^- e arctan. 2 ^-( 1 + / ^ 1 + / ^ 



1-7 2v^ 



1 



2/x(l + /) 



With our data, we get E = 0.99 ± 0.05 (la). The un- 
certainty in E is dominated by the errors in the galaxy 
position. If the mass distribution has the same ellipticity 
as the light (e = 0.33), an external shear of 7 min = 0.21 
is needed to keep E inside the la bounds. Because of the 
absence of very close galaxies or a rich cluster in the field, 
we do not expect such a large shear. For an ellipticity 
of e = 0.2, the minimal shear decreases to 7 m ; n = 0.11. 
Fig. [I| shows the possible parameters consistent with the 
measured positions and the flux ratio. Even for different 
parity of the images some very symmetrical models lead 




Fig. 4. Model parameters e and 7 consistent with the ob- 
servations. Models with more than two images are shown 
hatched. The dotted lines are la bounds of the measure- 
ments, the dashed line is a lower limit for 7 independent 
of f x . The solid line shows a symmetrical configuration 
with f x = 1. 

to more than two images (up to eight images are possible 
for SIEMD+shear models). 

Only rough estimates for the absolute amplifications 
can be determined from the observations. For a best-fit 
model, we get Ma « 70. Considering the errors, a lower 
limit for \M^\ of 27 (68% confidence) can be obtained. 

To estimate the mass and velocity dispersion of the 
galaxy, we use a spherical model (e = 0). For lens redshifts 
of z d — 0.3 (0.5), the mass inside the Einstein radius is 
M = 1.5 (2.4) x 10 11 h§Q Mq , and the velocity dispersion 
a v = 205 (230) kms" 1 (Q = 1, A = 0). The implied mass- 
to-light ratio is M/Lr = 5 (2) in solar units. 

The expected order of magnitude for the time delay is 
about weeks. A better estimate must await more stringent 
constraints on the galaxy position. Given the geometry of 
the system, off-center spectroscopy of the galaxy should 
be possible from the ground under excellent seeing condi- 
tions, or with STIS onboard the HST. 
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